Mean-field magnetization relaxation in conducting ferromagnets 
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Collective ferromagnetic motion in a conducting medium is damped by the transfer of the mag- 
netic moment and energy to the itinerant carriers. We present a calculation of the corresponding 
magnetization relaxation as a linear-response problem for the carrier dynamics in the effective ex- 
change field of the ferromagnet. In electron systems with little intrinsic spin-orbit interaction, a 
uniform magnetization motion can be formally eliminated by going into the rotating frame of refer- 
ence for the spin dynamics. The ferromagnetic damping in this case grows linearly with the spin-flip 
rate when the latter is smaller than the exchange field and is inversely proportional to the spin-flip 
rate in the opposite limit. These two regimes are analogous to the "spin-pumping" and the "breath- 
ing Fermi-surface" damping mechanisms, respectively. In diluted ferromagnetic semiconductors, the 
hole-mediated magnetization can be efficiently relaxed to the itinerant-carrier degrees of freedom 
due to the strong spin-orbit interaction in the valence bands. 



PACS numbers: 76.50.-|-g,75.45.-|-j,85.75.-d 

Relaxation of the ferromagnetic magnetization dynam- 
ics is well understood phenomenologically, being often 
accounted for by a single dimcnsionlcss parameter, the 
so-called Gilbert damping a} The equation of motion 
conserving the magnetization magnitude is written for 
the local magnetization-direction unit vector m as 



(?tm = —7m X Hcff + am x 9tm , 



(1) 



where 7 is the (minus) gyromagnetic ratio. The first 
term on the right-hand side describes the motion of 
the magnetization, M = Mgrn, in the effective field 
Hcff = — c^M^-iM], which preserves the magnetic energy 
i?[M] in the presence of applied, crystal, exchange, and 
demagnetization fields.'S. The second term characterizes 
the dissipation of the magnetic energy due to coupling 
with other degrees of freedom. In the case of small- 
angle motion near an equilibrium rotational-symmetry 
axis, Eq. describes a damped circular precession with 
frequency uj = 7i?efT- In the presence of anisotropics, 
both CO and a become tensor quantities, and the trajec- 
tories elliptic. For the purpose of our discussion, it is 
sufficient to treat the simple case of circular precession 
with a scalar damping a. 

Despite decades of experimental and theoretical stud- 
ies of itinerant ferromagnetism in metals and, more re- 
cently, in semiconductors, the microscopic origin of a is 
still not fully understood. One possible proposed mecha- 
nism involves a transfer of the angular momentum (and 
energy) of a nonequilibrium ferromagnetic configuration 
to the itinerant electrons via the exchange interaction, 
with a subsequent spin-orbit relaxation to the lattice. 
Such a process has been studied extensively within the 
s — d model, see, e.g., Refs. 00, although its imphed 
applicability to the itinerant transition-metal ferromag- 
netism has not been demonstrated. The s — d picture was 
resurrected recently- to address the question of magne- 
tization relaxation in the ferromagnetic semiconductor 
(Ga,Mn)As, where the ferromagnetism originates in the 
hole-mediated exchange interaction between the substi- 



tutional (paramagnetic) spin-5/2 Mn atoms ifi This letter 
puts forward a description of the magnetization damp- 
ing due to the exchange interaction between the localized 
magnetic orbitals and the itinerant carriers in ferromag- 
netic metals and semiconductors, by reducing the prob- 
lem to treating the carrier dynamics in a time-dependent 
uniform exchange field. 

There is at least one Gilbert-damping mechanism 
which can be identified and measured separately from the 
others in metallic ultrathin films. Its origin is nonlocal: 
The motion of a small ferromagnet with a large surface- 
to-volume ratio pumps spins into adjacent conductors, 
which can then be damped by spin-flip scattering out- 
side of the ferromagnet, leading to the Gilbert form of 
relaxation with an a which can dominate the intrinsic 
damping^. This process was thoroughly studied exper- 
imentally in ultrathin &lma&^ giving remarkable agree- 
ment with the parameters obtained by first-principles 
band-structure calculations.^*^ Bolstered by the success 
in understanding the enhanced damping of thin films 
by using a picture of the ferromagnetic relaxation via 
time-dependent exchange interaction with itinerant car- 
riers, we apply related ideas to formulate a framework for 
studying intrinsic relaxation of conducting ferromagnets. 

Consider an sp ~ d model of a conducting ferromag- 
net, where the spins J of the itinerant s or p orbitals 
(either electrons or holes) are polarized by an exchange 
field fl along the magnetization direction m of localized 
d orbitals^ii 



H{t) = n- nin{t) ■ 3 . 



(2) 



Here is a (time-independent) Hamiltonian which de- 
pends on the host band structure. The exchange field 
can be induced by the localized paramagnetic impuri- 
ties, such as substitutional Mn atoms in (Ga,Mn)As. Al- 
though such exchange can be highly nonuniform on the 
atomic scales, we are making a simplifying assumption in 
Eq. |2Jl of a uniform field whose magnetization direction 
m is treated classically in the mean-field approximation. 
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In particular, the magnetization is taken to be spatially 
uniform on the relevant length scales of the carrier dy- 
namics. We choose as a concrete example for some of our 
discussions the spherical Luttinger Hamiltonian for the 
spin-3/2 holes in the valence bands of a dilute p-doped 
semiconductor (e.g., GaAs, Si, or Ge): 

H-(2m,)-i{[7i + (5/2)72]p2- 272(p-JA)'} , (3) 

where rue is the free-electron mass and the 7^ are the so- 
called Luttinger parametersii^ The spin-orbit term cou- 
ples the hole momentum p with its spin J. For the valid- 
ity of the four-band model ((H)), the carrier density must 
be low enough that the Fermi energy is smaller than the 
intrinsic host spin-orbit interaction energy. For the dis- 
cussion of spin-1/2 electron systems, we set 72 = 0. Sup- 
pose the magnetization of the localized orbitals varies 
slowly in time (being uniform at all times), so that the 
time-dependent m modulates the Hamiltonian ||2J) adi- 
abatically. This means that the system equilibrates on 
time scales faster than the motion of m and all the quan- 
tities parameterizing the carrier Hamiltonian stay con- 
stant. Many-magnon scattering^^ is disregarded. Such 
a time-dependent long-range ferromagnetic order can be 
achieved in ferromagnetic resonance (FMR) experiments 
on thin films of transition metals^^ and semiconductors,^^ 
Consider the small-angle dynamics of the unit vector 
m(t) in Eq. |(2Jl near the z axis. Suppose the equilibrium 
value of the average spin is coUinear with the magne- 
tization. The variation Sm{t) = ni(i) — z will induce 
the uniform spin density 5jx{(jj) = Xj:cjj:i^)^^'>^xi^) + 
Xj^jy{uj)^lSmy{uj) along the x axis, and similarly along 
the y axis, where Xjj is the spin-spin response function 
at a finite fl (that is we are developing a perturbation the- 
ory for small variations in the direction of the exchange 
field, not assuming the smallness of its magnitude SI). 
For a system which is spin-rotationally invariant around 
the z axis, Xj^j^ = Xj j i which will be assumed in the 
following. The spin density j, in turn, corresponds to the 
effective field HcS = {^/Ms)j which gives a contribution 
to the magnetization equation of motion In the low- 
frequency limit, it translates into the damping coefficient 

a = {jiiyAQ lim Imx,^,. {u;)/lj . (4) 

The susceptibility Xj^j renormalizes the local-spin gyro- 
magnetic ratio 7 if the carrier-spin density is comparable 
to the local-spin density. Eq. I^J can also be obtained by 
equating the energy dissipated into the itinerant degrees 
of freedom by the moving magnetization and the work 
done by an rf magnetic field applied against the viscous 
Gilbert term in Eq. |p[J, at a steady magnetic preces- 
sion. Eq. Q is the most basic equation in this paper 
and may be taken as the definition of a. In the follow- 
ing, we formally evaluate a for electron and hole systems 
and discuss its dependence on the disorder composition. 

In the absence of spin-orbit interaction in the band 
structure, 72 = 0, the average spin density moves in the 




FIG. 1: Gilbert damping, Eq. in units of amax = 

7jo/(2Afa) as a function of the normalized spin-flip rate. In- 
set: Geometry of the model. 

exchange field fl as 

dtj = njx m(i) - [j - jom(t)] /n , (5) 

where the last term is a phenomenological relaxation due 
to impurities, characterized by the transverse spin-flip 
time T2. We assume here that m(t) undergoes a slow 
motion (on the scales set by frequencies fl and T^^)- It 
is convenient to transform Eq. 10) into the frame of ref- 
erence (for spin variables) moving together with m(i): 
If, for example, m is instantaneously rotating with fre- 
quency Lu around the y axis in the laboratory frame, see 
Fig. n inset, it is stationary in the rotating frame and 
there is a new (Larmor) term cijj x y on the right-hand 
side of Eq. which tries to polarize spins along the y 
axis. Because the motion is slow, we can solve for j as 
the (instantaneous) stationary state in the moving frame 
of reference. We then find the torque that j exerts on m 
to get 

« = -IXi^jy^/Ms = [1 + (r!T2)2] 7J0/M, , (6) 

where Xixiy stationary (real-valued) response func- 

tion in the rotating frame, for m pointing along the z 
axis. We have thus reduced the calculation of the time- 
dependent response in the laboratory frame, Eq. to 
the static response in the rotating frame. Such a trans- 
formation can be done in general for a spin-rotationally- 
invariant Hamiltonian with no spin-orbit interaction in 
the band structure. We plot Eq. ^ in Fig. ^ The equi- 
librium spin density jo can be calculated from the specific 
form of the Hamiltonian. a vanishes at both small and 
large spin-flip rates. 

The damping parameter © scales differently with the 
spin-flip rate depending on how it compares with the ex- 
change energy. The low spin-flip rate regime, a oc Tj"^, is 
completely analogous to the "spin-pumping" damping^ 
of thin fllms in contact with a "spin-sink" conductor: 
The moving magnetization pumps spins into the itiner- 
ant carriers at a constant rate, which are then relaxed 
with a probability oc before exchanging spins with 
the ferromagnet. The difference is that now the spins are 
"pumped" into the ferromagnet's own delocalized states. 
The other limit, a oc T2, is simple to understand since 
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~ jo{t) — '72i9tjo(i), in the laboratory frame when the 
relaxation rate dominates the dynamics of j (t) . a oc T2 
then follows from the torque oc (j — jo) x m. This is 
analogous to the "breathing Fermi-surface" mechanismi^ 
of the itinerant carriers, which try to accommodate the 
changing magnetization direction but are lagging behind 
with a delay time of T2 . In the presence of an anisotropic 
spin-orbit interaction in the crystal field of a metallic 
ferromagnet, the breathing Fermi surface gives an ad- 
ditional contribution to damping, which scales linearly 
with the band-structure relaxation timei^^ 

It is interesting to note that Eq. ® reduces to 
the random-phase approximation result for the long- 
wavelength magnon lifetime due to the s — d interaction 
with spin-1/2 conduction electrons, which was obtained 
in Ref. fusing a fully quantum- mechanical description: 

a = [1 + {nT2f] 7 [nm*kF/{^n'^h)\ /Ms , (7) 

where m* is the band-structure mass and is the Fermi 
wave vector, and it was assumed that Ml <C Ep (the 
Fermi energy). The quantity in the second square brack- 
ets is just the total carrier spin density. Eq. was 
used in Ref. to explain the measured damping in thin 
permalloy films, which scaled linearly with the film resis- 
tivity, as expected due to the Tj"^ prefactor in the rele- 
vant limit of a large exchange energy, Vl ^ T^j"^' i'^ the 
transition-metal ferromagnets. Unlike the case of the fer- 
romagnetic semiconductors, the direct application of the 
s — d model result is however questionable in the case of 
itinerant ferromagnetism in the transition metals where 
the separation between the magnetic and conducting or- 
bitals is unphysical. 

We now turn to a discussion of the application of these 
results to magnetization relaxation in hole-doped mag- 
netic semiconductor (Ga,Mn)As. Let us first make a 
rough estimate of the damping coefficient using Eq. ((HJ: 
The largest achievable value of amax = 7jo/(2A^s) occurs 
when the holes are fully polarized giving amax ~ 0.1 — 0.3, 
roughly one third the ratio of the hole to the substitu- 
tional Mn concentrations. For realistic samples with a 
spin polarization of the order of unity, therefore, amax ~ 
0.1. The damping a is further suppressed by the factor 
a/amax = 2nT2[l + (f7T2)2]-i < 1. For clean bulk sam- 
ples of GaAs, the spin-flip relaxation time is ^ 100 fsiii 
For approximately 5 % Mn doping, hVt 0.1 eVji^ so 
that ^T2 ^ 10, puts one in the a oc Tj^^ regime with 
a ~ 0.01. Shorter spin-flip times would thus result in 
larger damping. Experimentally, the impurity scatter- 
ing is likely to be the easiest parameter to vary in order 
to "engineer" a desired a. For a bulk sample, the strong 
spin-orbit coupling 72, however, makes the validity of the 
phenomenological equation |SJ), and thus result lO, ques- 
tionable. Besides, the crystal anisotropy would require a 
further refinement of the analysis. We thus have to return 
to Eq. Q in order to derive a reliable result. Evaluating 



the response function for a noninteracting Hamiltonian 
yields 

a = ^lini^^|(^|J,|j)|2[/(e,)-/(e,)] 

x6{huj + e, ~ ej) , (8) 

where /(e) is the Fermi-Dirac distribution and i, j label 
one-particle eigenstates in the sample of volume V . If the 
lattice vector k is conserved, ^j^- /V = / d^k/ (27r)'^, 
where a, h label spin states. For a perfect crystal, there- 
fore, a vanishes, as expected (unless there is a finite- 
measure Fermi surface area with a spin degeneracy) . The 
role of the relaxation on lattice defects was formally intro- 
duced in Ref. by broadened one-particle spectral func- 
tions, Aa(e) = r/[(e - eka)^ + r2/4], as follows: 

ct = -TT- lim - / 7 — > {ka\JJkb)r / - — -,r 
Ms ^^ou;J 2^ 3Z^'^ ' ' J 2^ 2 

ab 

X Aka(e)Akb(e -I- hiu) [/(e) - /(e + huj)] , (9) 

which was obtained by evaluating the local-spin suscep- 
tibility after integrating out the itinerant-carrier degrees 
of freedom. They^ find a nonmonotonic behavior of a 
as a function of the phenomenological scattering rate F 
for a realistic (Ga,Mn)As band structure: a oc F"^ as 
F — > (after taking the uj limit first) and, after 
passing through a minimum, a increases monotonically 
with large F. Eq. ^ however appears to have a problem 
for the large momentum-scattering rate, F, asymptotic: 
When F ^ 72, the D'yakonov-PereliS spin-relaxation 
rate for Hamiltonian Q scales as 7|/r; the spin-spin 
response corresponding to Eq. (O, on the other hand, 
has a 1/F cutoff in the time domain, resulting in a spin- 
relaxation rate growing linearly with F. 

In our discussion, we have assumed that the ferromag- 
netic magnetization is moving without specifying the ex- 
act mechanism of how the motion is initiated. Doing 
this by, e.g., applying an external magnetic field (with a 
large dc and small rf components) will of course affect 
the form of the Hamiltonian (0) for the itinerant carri- 
ers. Our results for the Gilbert damping will stay un- 
affected, however, as long as the exchange energy HQ is 
much larger than the carrier Zeeman splitting in the ap- 
plied field, and the ferromagnetic magnetization is mostly 
supplied by the localized orbitals (otherwise one has to 
take into account the energy pumped by the rf field into 
the carrier-magnetization dynamics) . 
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